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On 2010 October 13, the X-ray astronomical satellite Rossi XTE, during the observation of the 
newly discovered accretion powered X-ray pulsar IGR J17480-2446, detected a lunar occultation of 
the source. From knowledge of lunar topography and Earth, Moon, and spacecraft ephemeris at the 
epoch of the event, we determined the source position with an accuracy of 40 mas (la c.l.), which 
is interesting, given the very poor imaging capabilities of RXTE (~ 1°). For the first time, using 
a non-imaging X-ray observatory, the position of an X-ray source with a subarcsecond accuracy is 
derived, demonstrating the neat capabilities of a technique that can be fruitfully applied to current 
and future X-ray missions. 

Subject headings: Moon - pulsars: general - pulsars: individual (IGR J17480-2446) - stars: neutron 
- X-rays: binaries 



1. INTRODUCTION 

Since the dawn of the X-ray astronomy, the lunar 
occultation technique has been used to study th e po- 
sition and structure of X-ray sources (see, e.g., iBornl 
119791 and references therei n). An X-ray mi ssion like 
the ESA mission EXOSAT (|Tavlor et al.lll98l was con- 
ceived to have, as part of the core progra m, the study 
of X-ray sources using lunar occultations (jBornl 119791 : 
iBorn fc Debruniierlll979h . However, the development of 
X-ray mirrors with arcsecond resolution, rendered, de 
facto, ineffective such possibility. Indeed, the precision 
on the spacecraft ephemeris and the lunar surface to- 
pography at that time, gave an uncertainty greater than 
the i maging resolution of such optics. A decade passed 
before [Mereghetti et al.l (|1990T ) proposed again this tech- 
nique for missions like XMM-Newton, but there is no 
trace in the last 20 years of literature of such an appli- 
cation. 

IGR J 17480-2446 was detected for the first time on 
2010 October 10 b y the IBIS-ISGRI instruments on 
board INTEGRAL (jBordas et al.l l20Toh in the globu- 
lar cluster Terzan 5. While observed by the X-ray 
observatory Rossi XTE (RXTE hereafter), the source 
showed an eclipse, first attributed to the companion star 
(|Strohmaver fe Markwardt 2010;). A further investiga- 
tion revealed that a very rare and serendipitous event 
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happened: the sour ce was eclipsed by the Moon while 
o bserved by RXTE (iStrohmaver et all [20Toh . 

Pool ev et al.l ()2010 ) . on the basis of a Chandra obser- 
vation of this outburst, gave the most accurate source 
position, identifying IGR J17480-2446 with a quiescent 
LMXB pre viously detected in a Chandra obs e rvatio n of 
Terzan 5 (|Heinke et all 12001) ■ iTesta et al.l (|2^TTI) re- 
port the identification of a possible near-IR counterpart 
through observations of the field with Adaptive Optics 
systems before and after the discovery of the source. A 
detailed discussion will be presented in a forthcoming pa- 
per (V. Testa et al . 2012 in preparation). Very recently, 
iPatruno fe Milonel (|2012f) report the identification of the 
candidate counterpart of IGR J 17480-2446 in the opti- 
cal bands in Hubble Space Telescope 2003 archival obser- 
vations whose positi on is in agreement with the X-ray 
positions reported bv lPoolev et al.l ([2010D and this work. 

In the following we will show how, taking advantage 
from this fortunate event, it is possible to determine the 
position of IGR J17480-2446 with a subarcsecond preci- 
sion. The discussion is divided in two parts, where we 
explain the occultation technique and provide a detailed 
analysis of the error budget, respectively. 

2. THE TECHNIQUE 

The precision of the lunar occultation technique de- 
pends on four factors: (1) the precision on the position 
of the Moon, (2) the precision on the spacecraft position, 
(3) the precision on ingress and egress epochs, (4) and 
the precision of the lunar surface topography. 

During the ingress and egress of a lunar eclipse, the 
occulted X-ray source lays on the projection of the lunar 
rim on the plane of sky as seen by the telescope (in this 
case Rossi XTE). The source position in the plane of sky 
is then given by the intersection of these two profiles. In 
general, the possible intersections are two, one of which 
can be usually excluded, under some conditions, if a prior 
estimate of the source position is available. 

To reconstruct the profile of the lunar rim in the plane 
of sky, as seen from RXTE at a given epoch, we first con- 
sidered the position of the spacecraft and of the Moon 
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R.A. - 17 h 48 m 04*8245 (deg) 

Figure 1. Moon topography projection on the plane of sky as 
seen by Rossi XTE at eclipse ingress (red curves) and egress (green 
curves). Each curve is the projection on the plane of sky of a lunar 
profile at constant selenographic longitude. The envelope of these 
curves is the lunar rim as seen by RXTE. In the inset the detail 
of the intersection between the two profiles d etermining our best 
source position is reported. Circles identify the Hcinkc ct al. (2006) 
Chandra position and the position derived in this work (black lo- 
ci., gray 2<r c.l.), respectively. 

barycenter at a given epoch (in our case ingress or egress 
epochs). The space-time reference frame we adopted is a 
Cartesian coordinates in a geocentric international celes- 
tial reference frame (ICRF/J2000.0) for the position and 
terrestrial time (TT/TDT) for time. 

Moon position (t) and spacecraft position txteW 
are available as astrometric positions, i.e., as a function of 
time in a reference frame centered on the Earth barycen- 
ter. Since the Moon astrometric and apparent positions 
(as seen from RXTE) differ because of the travel time of 
X-ray photons to cover the Moon- RXTE distance d, we 
must correct for this time, At = d/c, where c is the speed 
of light. The expression to evaluate the vector distance 
d is reported below 

d = r XTE (i) - [r<t (t - At) - (r e (i) - r®(* - At))] , 

w 

where Tc^,(t) is the Earth position with respect to the 
solar system barycenter. The expression between square 
bracket represents the Moon position at time t — At in 
the reference frame at time t. Since At is a function 
of d, Equation([T]) has to be solved iteratively to find d. 
We used a fixed point method which converged to the 
required accuracy within few steps. 

Since all the available lunar surface topographic maps 
are expressed with respect to the selenographic reference 
frame, the next step was to place the selenographic co- 
ordinate system in ICRF/J2000.0 reference frame. We 
were then able to project the lunar surface topography 
as seen by RXTE on the plane of sky at both eclipse 
ingress and egress epoch (see Figure [T]). The coordinates 
of the intersection point (inset in Figure [l} are reported 
in TableHJ 

3. OBSERVATIONS AND DATA ANALYSIS 



The source was discovered on 2010 October 10 
(jBordas et al.ll2010D . and RXTE started the observa- 
tion of the source on October 13. The lunar eclipse is 
present in the first pointed observation of IGR J17480- 
2446 (Obsld 95437-01-01-00) performed by the Propor- 
tional Counter Array (PCA)/ RXTE instrument. The 
data are available in good-xenon packing mode, with 
maximum temporal (1/is) and energy resolution (256 
channels). For the kind of analysis we performed, pho- 
ton arrival times in data where not reported to the solar 
system barycenter, since we are interested in the instant 
at which photons reach the spacecraft. Fine clock correc- 
tions were not applied since they have negligible effect in 
the present analysis. Indeed, their magnitude (< 60/is) 
is two orders of magnitude smaller than the uncertainties 
on ingress and egress epochs we obtained. 

In the following, a detailed discussion of all the uncer- 
tainties is given. 

3.1. The Moon 

In the determination of the source position, the sources 
of error related to the Moon are two: the positional er- 
ror of the Moon barycenter and the uncertainty on the 
knowledge of the Lunar topography. 

3.1.1. Uncertainty on the Moon Position 

In our work, we adopted the Moon position with 
respect to the Earth as reported in DE421/LE 421 
ephemeris (|Folkner et al.ll2008t iWilliams et al.ll2008f ). 

The uncertainty at la confidence level (hereafter c.l.) 
in the lunar orbit for DE421 (taking into account possible 
systematic errors, not just formal uncertainties from the 
least-squares fit) is about 5 m in right ascension, 2.4 m 
in declination, and 0.5 m in range (W. F. Folkner 2012, 
private communication). 

The position of Earth and Moon bary centers, and the 
spatial orientation of their spin axis as a function of time 
was obtained using the JPL's NAIF SPICE Toolkit We 
sampled these positions every 15 s. 

During the iterative process described in Section [3J we 
interpolated between these positions with a cubic spline. 

3.1.2. Uncertainty on the Moon Topography 

The topography of the lunar surface was measured by 
Earth-observation and lunar missions. The most accu- 
rate topographic maps of the moon was obtained by laser 
ranging altimetry made by probes orbiting the Moon. 

The very first one of such a topographic lunar map was 
done by the NASA Clementine mission (|Nozette et al.l 
19941) . with a radial precision of ~ 40 m. In this work, 
we analyzed the topographic data from the laser altime- 
try instrument LALT on b oard the Japanese lunar ex- 
plorer Kaguya (SELENE) (lAraki et al.l 120091 [20Toh . in 
particular the global grid topographic data of the Moon 
(LALT_GGD_NUM, ver. 1, hereafter LALT data). The 
LALT data are referenced to the sphere of 1737.4 km ra- 
dius based on the gravity center of the Mean Earth/Polar 
Axis body-fixed coordinates of the Moon0 The LALT 
data grid resolution is 0.0625 (1/16) deg (|Araki et al.l 

8 Version N0064 http:// naif . jpl .nasa. gov/naif /toolkit .html 

9 For details on the lunar coo rdinate system adopted, see 

http : //lunar . gsfc .nasa. gov/library/LunCoordWhitePaper- 10-08 .pdf 
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Figure 2. Light curve of PCA data around the eclipse of 
IGR J17480-2446. The light curve is binned at 1 s. In the two 
insets the light curve of eclipse ingress (left) and egress (right) 
binned at 10 ms are reported. 



[20091 [20101 iFok et alJl20lTl). The radi al topographic er- 
ror is 4.1 m (la c.L. IAraki et al.ll2009D . 

3.2. Uncertainties on Ingress and Egress Epochs 

During the observation in which the lunar occultation 
occurred, the source flux was nearly constant with the 
exception of a type I burst (outside the time interval 
shown in Figure [2]). Excluding the type-I burst, the av- 
erage count rates outside and during the eclipse where 
254.3(9) and 36.2(9) counts s _1 , respectively (number in 
parenthesis indicates ltr errors on the last digit). These 
values are obtained fitting the light curve (binned at 1 s) 
inside and outside the eclipse with a constant models. 

To obtain an accurate determination of the epoch of 
eclipse ingress and egress, we binned the light curve at 
100/is (on such short time intervals no more than one 
event occurred in each bin, implying that the events can 
be considered unbinned, in practice) and fitted separately 
eclipse ingress and egress with step (Heaviside) functions. 

Due to the photon paucity, a least-square method was 
not suitable to fit the ingress and egress epochs (see insets 
in Figure [5]) . We therefore used an unbinned maximum 
likelihood method, which is m ore suitable in such cases 
(|Bevington fc Robin son 2003) . We have taken advantage 
from prior knowledge of the count rate inside and out- 
side eclipse, so that the fit is mono-parametric with the 
only parameter the transition time. We fitted separately 
ingress and egress, considering for each fit the data in 
a 10 s time interval centered on ingress (or egress). To 
estimate the uncertainty to be associated to the method, 
we simulated 10,000 data set with the same characteris- 
tics of the original one, and studied the distribution of the 
best-fit values. The confidence interval obtained is asym- 
metric, reflecting the large difference between statistics 
inside and outside eclipse. The la uncertainty results to 
be 7 ms on the high count rate side, while is 9 ms on the 
low count rate side. The ingress and egress epochs and 
the associated la c.l. are reported in Table [TJ 

3.3. Uncertainties on Spacecraft Position 

A complex part of the analysis was the determina- 
tion of the Rossi XTE spacecraft positional errors. In 
the following, a brief description of how the spacecraft 
ephemeris are derived is given. 

The Rossi XTE spacecraft positions (and velocities) 



Table 1 

Position of IGR J17480-2446 and Times of Eclipse 
Ingress/Egress. 



Parameter 



R.A. (ICRF/J2000.0) 17 h 48 m 04?8245(26) 
Decl. (ICRF/J2000.0) -24° 46' 48'.' 88(4) 



T in (MJD TT) 



55482.03581125 



+ 10 



T cg (MJD TT) 55482.04173510^0 



Note. — Numbers in parenthesis are la c.l. errors on the last 
digit(s). 



are stored in the orbit files as Cartesian co ordinates with 
a sam pling rate of 60 s. As reported in Uahoda et al.l 
(2006), the position is obtained from the spacecraft 
ephemeris, estimated by the Goddard Flight Dynamics 
Facility by fitting the spacecraft position on a 2/3 days 
basis. This procedures implies that orbit ephemeris of 
contiguous orbitfiles are more or less correlated. In each 
orbit file, the spacecraft best-fit ephemeris is also extrap- 
olated forward by 10 hr, overlapping with the next day 
solution. The overlap between two contiguous orbit files 
is not perfect, reflecting the orbit variations due to the 
solar radiation pressure, and the statistical uncertainties 
in the determination of the spacecraft ephemeris. To 
estimate the stat istical uncerta i nty on the spacecraft po- 
sition, following Uahoda et al.l (|2006[ ) , we analyzed the 
distribution of the distance between the spacecraft posi- 
tions reported in the overlapping time intervals from the 
orbit number 5000 (54354 MJD) to orbit number 6203 
(55556.0 MJD), covering a time span of 3.2 yr around 
the eclipse epoch. This distribution shows three peaks, 
centered on 0, ~ 0.01 m, and ~ 50 m, respectively, re- 
flecting the aforementioned degree of correlation between 
contiguous orbit files. 

In line with what is discussed in Uahoda et al.l (2006) 
and following a discussion with C. B. Markwardt (2012, 
private communication), as an upper limit on the un- 
certainty on the spacecraft position, we considered the 
statistical distribution of the ~ 50 m peak. 

From this distribution we derive the following confi- 
dence intervals for the Rossi XTE position in the con- 
sidered time interval: 48 m at la, 121 m at 2a, 311 m at 
99%, and 970 m at 3a. 

4. DISCUSSION AND CONCLUSIONS 

As described in Section^ we were able to determine 
the position of the source as intersection of the lunar rim 
as seen by RXTE at eclipse ingress and egress (see Figure 
[T|). In the following, to calculate the overall uncertainty 
of IGR J17480-2446 position, we consider all the source 
of error as independent. The sources of uncertainty can 
be subdivided in two kind: temporal (ingress and egress 
epoch, fine clock corrections) and spatial (spacecraft and 
lunar position, lunar topography). 

The position on the plane of sky of the source is given 
by the direction of the segment joining the spacecraft 
and the point on the lunar surface in which the eclipse 
ingress and egress of the X-ray source occurred. This im- 
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Table 2 

Error Budget 



Parameter 


A H.A. (mas) 


A Decl. (mas) 


Moon cphcmcris 


±2.7 


±1.3 


Rossi XTE ephemeris 


±26 


±26 


Eclipse ingress and egress epochs 


+29 
-25 


+24 
-21 


Moon topography radial precision 


±2.7 


±2.7 


Grand total 


+ay 

-36 


+3b 
-33 



Note. — Statistical errors are intended to be at lcr confidence 
level. 

plies that the angular uncertainty on the direction of such 
a segment is obtained from the ratio of the overall posi- 
tional uncertainty on a plane perpendicular to the source 
direction and the Moon-spacecraft distance. The overall 
uncertainty on the plane of sky is obtained by summing 
in quadrature all the spatial uncertainties involved. 

As reported in Section 13.11 the positional uncertainty 
on the position of the Moon barycenter are 5 m in 
R.A., and 2.4 m in decl. The corresponding angular 
uncertainty, calculated considering an average Moon- 
spacecraft distance of 3.844 x 10 5 kmJ3 is of 2.7 mil- 
liarcseconds (hereafter mas) in R.A. and 1.3 mas in dec- 
lination. The error associated to the lunar topography 
adopted is the radial uncertainty of 5 m. The angular 
uncertainty we obtain is 2.7 mas in both axes. 

Using a Monte Carlo technique, we estimated a confi- 
dence interval of ±7 ms and l~g ms for ingress and egress, 
respectively. To convert these temporal uncertainties 
in spatial uncertainties on a plane perpendicular to the 
source direction, we multiply them by Moon-spacecraft 
relative speed projected on the plane of sky v re i. The 



„R.A. _ 



5.70 kms" 



„Deci. 



= 4.71 



moduli of these velocity vectors are: 
vP ecL = 2.26 km s~ 

km s _1 , where the subscripts in and eg indicate ingress 
or egress epoch, respectively, while the superscripts RA 
and Dec indicate the corresponding axis in the plane 
of sky. Proceeding as described above, we obtain for 
R.A. and decl. the following lcr confidence intervals: 
Aa in = ±21 mas ; Acicg = +\l mas, A(3 in = tn mas, 
Aficg = mas. The uncertainties in ingress and egress 
where summed in quadrature. 

Finally we considered the uncertainty in the spacecraft 
position, which is 48 m at lcr confidence level, corre- 
sponding to an uncertainty of 26 mas on both axesF 1 ! 

To summarize, we report in Table [2] all the uncertain- 
ties discussed above at lcr c.l. on R.A. and deck, together 
with the grand total, obtained summing them in quadra- 
ture. 

Our result con firms the association made by 
iPoolev et al.l (120101) of a quiescent X-ray source observed 



Heinke et al. H20fJ6l) in the globular cluster Terzan 



by 

5. Indeed our posi tion is O'.'l from th e quiescent X-ray 
source detected bv lHeinke et al.l (|2006f ) in a 40 ks Chan- 
dra X-ray observation of the globular cluster Terzan 5. 

10 The Moon— spacecraft distance is 383418.495 km at ingress, 
and 385476.095 km at egress. In our calculations, we considered 
the average. 

11 We note that, formally, we should project the spacecraft posi- 
tion uncertainty on the plane of sky. Such a procedure implies that 
the uncertainty to be associated should be a factor smaller. 
However, to be conservative, we prefer to use the original uncer- 
tainty. 



No te that the accuracy on the source position reported 
bv iHeinke etall (|2006r ) is 0"06 (lcr c.l.), slightly worse 
than the uncertainty derived form the Moon occultation 
(~ 0'.'04). Moreover, V. Testa et al. (2012, in prepara- 
tion) use the high precision on the X-ray source position 
obtained with Moon occultation to identify a possible 
near-IR counterpart even in a extreme crowded field of 
view as the center of a globular cluster. 

In the following we discuss the applicability of this 
method to present and future X-ray missions. In par- 
ticular we briefly compare RXTE with Chandra, XMM- 
Newton, NuSTAR, and LOFT discussing the positional 
accuracy obtainable with this technique. As discussed in 
Section[2j the accuracy on the source position depends on 
the accuracy in the ingress and egress epochs, which in 
turn depends on the source and background count rates 
through the following relation St ~ (1 + i?bkg/ -Rsrc) /-Rsrc 
(at lcr confidence level), R src and i?bkg are the count rate 
of the source and the background, respectively. There- 
fore the total uncertainty, at lcr level, in source position, 
ag, can be written as 



d 2 



1 + -Rbkg/ 
Rsrc 



spc 



d? 



(2) 



where and cr spc are the uncertainties in the Moon 
rim and spacecraft position, d is the Moon-spacecraft 
distance. Since solar system ephemeris, and in partic- 
ular Lunar ephemeris, are already at an a ccuracy of 5 
m, w hile present lunar mission LRO/LOLA ()Smith et al.l 
I2010T ) is mapping its surface with an overall accuracy < 1 
m, the precision is mainly limited by the uncertainties on 
spacecraft position and on the ingress and egress epochs. 
The first one can easily be improved in present and fu- 
ture spacecrafts equipped with GPS. This means that 
the leading term in the uncertainty ag is the uncertainty 
on the ingress and egress epochs. For source dominated 
targets (i?bk g /-Rsrc << 1), the positional accuracy ob- 
tainable with this technique scales as the inverse of the 
source count rate, which, for a given spectral shape, de- 
pends on the effective area and energy band of the X-ray 
satellite. The effective area at 5 keV for present X-ray 
satellite ar e 1000 cm 2 for each Proportional Counter Unit 
in RXTE (jBradt et al][l99l) 800 cm 2 for EPIC/pn in 
XMM-New ton dJansen et al.ll200lD . 300 cm 2 for Chan- 
dra /EPIC (jWeisskopf et al.1120021), while the effectiv e are 
of NuSTAR will be 700 cm 2 (jHarrison et all [20ToT) an d 
of LOFT/LAD will be 110,000 cm 2 (jFeroci et al.|[20ll . 
Unfortunately, such a technique cannot be used with 
XMM-Newton, since it has a constraint on the angle 
between the pointing direction and the Moon (> 22°), 
because of the limitations imposed by the Optical Mon- 
itor. For bright sources, continuous clocking mode must 
be used with Chandra, in order to avoid pile-up photon 
losses. A future X-r ay mission with a ~10 m 2 effective 
area such as LOFT (|Feroci et al.l 1201 lT ) will then permit 
to go at a level below 5 mas. 

Note added in proof We have re-calculated the position 
of the 11-Hz pulsar ICR J 17480-2446 in Terzan 5 using 
data from LRO /LOLA for the Lunar topography. These 
data have an intrinsic uncertainty of 1 m on the Lunar 
surface. In this case we find the following source posi- 
tion: 17 h 48 m 04^8216(26), -24° 46' 48"88(4). This posi- 
tion differs from the one obtained from Kaguya/LALT 
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Figure 3. Moon topography of the 8° X 8° regions around the eclipse ingress (left panel) and egress (right panel) points. Each map is 
centered on the direction of the moon radius orthogonal to the source direction. The white circles identify the points on the lunar surface 
where the ingress and egress took place. The contour interval is 500 m. The white crosses indicated the most important lunar craters in 
the regions. 



by 43 mas, which is about 1 sigma the error box. 
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